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We analytically and numerically characterize the structure of hard-sphere fluids in order to review various 
geometrical frustration scenarios of the glass transition. We find generalized polytetrahedral order to be 
correlated with increasing fluid packing fraction, but to become increasingly irrelevant with increasing dimen- 
sion. We also find the growth in structural correlations to be modest in the dynamical regime accessible to 
computer simulations. 
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Two key problems in glass physics can be formulated 
as follows: (i) What is the mechanism that thwarts crys- 
tallization of a liquid and allows for the development of 
glassy behavior? (ii) Can one assign the fast increase of 
the relaxation time and the viscosity as one approaches 
the glass transition to a collective phenomenon and the 
growth of static correlations? Possible answers to both 
questions have been proposed, based on the concept of 
"frustration" . Frustration in this context refers to an in- 
compatibility between some local order that is preferred 
in the fluid and a global periodic tiling of spaced— . 

A detailed description of frustration requires knowl- 
edge of the relevant type of local order, a difficult (and 
unsolved) task in the case of a generic molecular glass- 
forming liquid. Hard-sphere fluids on the other hand are 
well suited for the task, because they are minimal model 
glass formers and their geometry is unambiguous. Frus- 
tration is then purely geometrical. For monodisperse and 
moderately polydisperse hard-sphere fluids, the locally 
preferred order is based on tetrahedral arrangements of 
4 particles, which when combined in larger units leads 
to icosahedral, or more rigorously polytetrahedral, order 
whose extension to the whole space is indeed frustrated. 
As often in physics, insight is gained by enlarging the 
parameter space. Geometrical frustration can be varied 
and fruitfully studied by varying the properties of the 
embedding physical space, either by introducing curva- 
ture^^ or by changing the number of dimensions. We 
follow here the latter route. 

Three-dimensional tetrahedral order generalizes to 
"simplicial order" in all dimensions, the simplex being 
the ideal packing of d + 1 particles in d dimensions: a 
triangle in 2 dimensions, a tetrahedron in 3 dimensions, 
a "hyper-tetrahedron" in higher dimensions. In d = 2, 
the simplicial order is not frustrated, but it is frustrated 
in d = 3, and even more so in higher dimensions. Di- 
mensional studies have already unambiguously answered 
the above question (i) in the case of hard sphere sys- 
tems&i, by showing that the local orders of the fluid and 



the crystal become increasingly different as dimension 
increases and that, as a result of this frustration, crys- 
tallization is strongly suppressed in dimensions d > 4. 
Whereas crystallization kinetically interferes in three- 
dimensional monodisperse hard-sphere fluids 8 , so that 
one needs to consider binary or polydisperse mixtures in 
order to study the metastable fluid branch^—, going to 
higher dimensions sidesteps this particular issue by sup- 
pressing crystallization even in monodisperse systems. 

The question that we address in this paper is then the 
second one, namely the connection between the dynam- 
ical slowdown and the growth of static structural cor- 
relations, which result in the super- Arrhenius, or frag- 
ile, dynamical behavior of glass-forming fluids. If indeed 
the phenomenon is collective, it should be accompanied 
by the development of nontrivial correlations to which 
one should be able to associate one or several typical 
length scales. It has been realized in the past two decades 
that slow relaxation in glass formers has an increasingly 
heterogeneous character, which can be characterized by 
an appropriately defined "dynamical" length obtained 
through multi-point space-time correlation functions 12 . 
The somewhat different issue that we consider here is 
whether the slowing down itself, as described by the in- 
crease of the structural relaxation time or the viscosity, 
is accompanied by the growth of a static length. 

There has been a number of proposals for relating dy- 
namics and structure that include some measure of the 
spatial extension of a characteristic locally preferred ar- 
rangement^— ii2r— . We investigate in the present work 
various measures of frustrated simplicial order in hard- 
sphere fluids in dimensions ranging from 3 to 8. Our 
results, which combine analytical and numerical ap- 
proaches, show that the growth of associated structural 
correlations under compression is very limited in the den- 
sity range that is accessible to computer simulation. This 
observation can be rationalized by invoking the strong 
degree of frustration that is present, even in d = 3, and 
that limits the putative collective behavior due to the ex- 
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tension of the frustrated local order. It could of course be 
that the structural correlations related to simplicial order 
do not capture the relevant static correlations that affect 
the dynamical slowdown. We therefore also investigate 
"order-agnostic" lengths, which characterize static point- 
to-set correlations; they are expected to provide upper 
bounds to the more specific length scales associated with 
local order while being more directly related to the co- 
operative nature of the dynamics^ 3 -. We also find that 
the evolution with density of such point-to-set lengths 
is quite limited in a domain where the relaxation time 
nonetheless increases by several orders of magnitude and 
the length associated with the dynamical heterogeneity 
grows by a factor of 4 or 5. 

The paper is organized as follows. In sections III and ITTl 
we review notions of higher-dimensional simplicial geom- 
etry as well as various geometrical proposals for linking 
the local geometrical order and the dynamical slowdown 
of fluids, respectively. In section IIII1 we obtain analyti- 
cal bounds for the fluid order and numerically character- 
ize the simplicial structure of hard-sphere fluids in d=3 
to 8. In section IIV1 we describe features of an "order- 
agnostic" structural length, calculate its growth for two 
hard-spheres fluids in c?=3, and compare the results with 
a bound extracted from the fluid dynamics. A brief con- 
clusion follows. 



I. GEOMETRICAL TILINGS OF SIMPLICES 

In order to describe ideal tilings of hard spheres we 
need to review some notions of Euclidean space R d , hy- 
perbolic space H d , and spherical space S d geometry. 
Note that the sphere space of (surface) dimension d, S d , 
is defined by the set of all vectors of equal length in 
R d+1 . Hence S 2 is the standard sphere in R 3 . Consider a 
(non-overlapping) configuration of N hard spheres in d- 
dimensional space of volume V— , at positions ri, . . . , tn. 
The neighbors of a particle are defined by the Delaunay 
tessellation associated with those positions, and the dis- 
tances are measured in units of hard-sphere diameter a. 
In R d , a configuration therefore has a number density 
p = N/V that is related to its packing fraction ip via 

<p = V d (a/2)p, (1) 

where V d (R) = R d -K d l 2 /T(d/2 + 1) is the volume of a 
d-dimensional ball of radius R. 

A simplex contains d+l vertices, and each of its faces is 
itself a lower-dimensional simplex (triangle in <i=2, tetra- 
hedron in d=3, . . .). It is thus the simplest polytope in 
any dimension. It also plays a central role in our analy- 
sis because, with probability one, simplices are the basic 
tiles of a Delaunay tessellation. Consider the i-simplices 
<7j(l), . . . , (Ji(iVj) of a Delaunay tessellation. Of course, 
Nq = N is the number of particles, N% is the number of 



bonds, and in general 

d 

YX-l) i N i = X , (2) 

where the Euler characteristic, depends on the topol- 
ogy of the space considered. Let 

qj,d(i) = #crd wrapped around <Jj(i), 

and 

1 Ni 

q ^ d = ]v"H*' rf W 

i=l 

be its average. We thus have that (70,2 is the number of 
triangles sharing a common point, and that go, 2 is the av- 
erage number of triangles per vertex in a configuration. 
In this work, we pay particular attention to bond spin- 
dles>&, q~d-2,d- We also consider Zd(j), the coordination 
number of the j th vertex, i.e., Zd(j) is the number of o\ 
wrapped around coG), and therefore 

1 N 

2d = jjH Zd ^- 

1=1 

In order to describe the regular tessellations of higher- 
dimensional space, we introduce Schlafli's notation. In 
this recursive description of generalized regular poly- 
hedra, the polytope {a\,...,a d } is composed of poly- 
topes {ai, . . . , ad-i} with vertex figure {02, . . . , ad} (see 
Ref. |26| . Chapter 7) . The vertex figure of a d-dimensional 
polytope is a [d— l)-dimensional polytope whose vertices 
are obtained by taking the middle of each edge emanat- 
ing from a vertex. The basic example of this inductive 
scheme is {p}, the regular polygon of p sides in R 2 . We 
then have the regular polytopes in R 3 , which are rep- 
resented by {p, q} for a polytope having faces made of 
{p}'s and vertex figures made of {gj's. Note that be- 
cause the polytope is regular, one can inscribe all the 
vertices on a sphere and thus {p, q} gives a tiling of the 
sphere S 2 . For instance, the cube {4,3} is composed 
of squares (the regular polygon {4}) and the vertex fig- 
ure is {3}, a triangle. The cubic tiling {4,4} in R 2 is 
composed of 4 squares surrounding each vertex. Note 
also that simplices are denoted {3 d-1 }. (In Schlafli sym- 
bols, exponentiation is understood to mean repetition, 
not multiplication.) Of particular relevance to this pa- 
per is the relation between Schlafli's notation and the 
spindle coordination. For {ai, . . . , ad} that is a regular 
tiling in R d or S d , one can indeed show that qd-14 — a d 
(Appendix [A| . 

Any polytope {ai, . . . , ad} in R d+1 can be seen as in- 
scribed in a sphere S whose center is the center of mass 
of the polytope. Inflating the polytope so as to make 
it round produces a regular tiling of S . For instance, 
the hypercube {4, 3 d-2 } in R d is also a regular tiling of 
S^ 1 . The usual case {4,3} therefore gives qd-2,d = 3 
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Simplicial tilings 


in 


curvature 


{3,6} 


R 2 





{3, n} with n > 7 


H 2 


— 4arccosh 2 (2 _1 esc -) 


600-cell {3, 3, 5} 


S 3 


tt 2 /25 ~ 0.39 


octahedron {3 d " 1 , 4} 


S d 


tt 2 /4 ~ 2.47 


{3,3,3,6} 


H 4 


4arccosh 2 ( 1 + v ^) ~ 4.505 



TABLE I. Summary of the regular simplicial tilings in d > 2. 
The curvature reported in units of a. 



on S 2 . The same argument illustrates that for the tiling 
{4,3 d - 1 } of S d one has q = 3. 

We stress that in the case of hard sphere fluids, we are 
interested in tilings by simplices. A detailed review of 
the possible regular simplicial tilings in various spaces is 
available in the Supplementary Material of Ref. 27. Here, 
we only report the curvature of the space that would 
embed these tilings (Table fl} , and highlight a few key 
observations. 

• A regular tiling of tetrahedra is found on the sphere 
that inscribes the remarkably large d — 4 regular 
polytope {3,3,5}. The inscribing sphere's curva- 
ture is much smaller than that of the generalized 
octahedron {3, 3, 4}, the only other polytope exclu- 
sively made of simplices. 

• By contrast, on the d-dimensional sphere with d > 
3, {3 d_1 , 4} is the only way to obtain a regular tiling 
of simplices, which can only be achieved on a sphere 
of much higher curvature. 

• The regular simplicial tiling {3,3,3,6} is found on 
a hyperbolic space of curvature larger (in absolute 
value) than the curvature of the sphere that in- 
scribes the tiling {3,3,3,4}, while the hyperbolic 
tiling {3,3, 6}'s infinite edge length disqualifies it 
as a possible reference for Euclidean packings of 
spheres. 

It is also interesting to note that the densest known 
crystal lattice in d=8, the root lattice Eg, is singularly 
dense, which arises from it having a remarkably high frac- 
tion of regular simplices in its midst 28 . Each of the sites 
in the Eg lattice is indeed surrounded by 2160 generalized 
octahedra and 17280 simplexes. 



II. FORMULATIONS OF GEOMETRICAL 
FRUSTRATION 

Various geometrical proposals have been made to ex- 
plain the growth of a static correlation length in glass- 
forming fluids. These approaches depend on the geo- 
metrical frustration that arises when optimally local and 
global orders are incompatible. All of these approaches 
thus also rely, in one form or another, on an analogy with 



hard-disk packings in R 2 . Thue long ago showed that the 
densest possible R 2 packing of disks is the triangular lat- 
tice 2 ^. Interestingly, both the densest closed-shell cluster 
- seven disks forming a hexagon - and the perfect sim- 
plex - three disks forming an equilateral triangle - tile 
space into this lattice. This system is thus understood 
having no geometrical frustration. In fact even the ther- 
modynamic transitions between the fluid and crystal are 
either continuous or very weakly first orde r 30 ' 31 . 

The behavior of hard disks in R 2 contrasts with that of 
disks on negatively curved space 5,32 , of mixtures of dif- 
ferent size disks and spheres ^ 11 ' 18 ' 33 ' 34 , or of spheres in 
higher-dimensional spaces 6,7 . In that sense, geometrical 
frustration between the fluid and the crystalline order is 
ubiquitous. By severely inhibiting crystallization, geo- 
metrical frustration extends the fluid branch to deeply 
supersaturated fluids, which gives access to the dynami- 
cally sluggish regime and glass formatio n 32 ! 35 ' 36 . 

We also wish to determine if the dynamical slowdown 
of these supersaturated fluids has a straightforward ge- 
ometrical origin. Can simple geometry provide both a 
measure of a growing static length and a mechanism 
for the slowdown? Higher dimensions provide an addi- 
tional constraint. Because the glassy phenomenology of 
hard-sphere fluids is remarkably robust to dimensional 
change s 35 ' 36 , we expect a geometrical proposal to be gen- 
eralizable to arbitrary dimensions. We use this dimen- 
sional requirement to critically review various proposals 
for the geometrical origin to the glass transition. 



A. Disclination Lines 

Regular simplices are the densest structure for packing 
hard spheres in any dimension. In dense fluids, where free 
volume is an important contribution to the fluid free en- 
ergy, one would therefore expect simplices to have a grow- 
ing structural representatio n 37 ' 38 . For Euclidean space in 
dimension d > 3, however, perfect simplices cannot tile 
space without defects. In d = 3, a mix of tetrahedra and 
octahedra is needed to assemble the face-centered cubic 
lattice; in d = 4, only octahedra are needed to generate 
the densest known lattice^, D4, and kissing cluster. 

Shortly after noticing that a perfect tiling of simplices 
with (J1.3 = 5 for all bonds is possible in d — 3 on the rel- 
atively gently curved S 3 space that embeds the {3,3,5} 
polytope (Table H|, it was suggested that d = 3 fluid con- 
figurations can be usefully described as frustrated with 
respect to that order. The fluid should then be under- 
stood as weakly defective with respect to that perfect 
tiling 1 ^. The structural defects that result from "flatten- 
ing" S 3 tilings to R 3 are known as disclinations. They 
are best understood by analogy to d = 2 tilings. Each 
particle in a perfect d = 2 triangular lattice of disks on R 2 
is part of six triangles, i.e., go, 2 = 6. Curving that plane 
results in irreducible disclinations that sit on disk cen- 
ters, and whose coordination differs from six. Irreducible 
disclinations with go, 2 > 6 are obtained in spaces of neg- 
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ative curvature^, e.g., H 2 , and with q Q , 2 < 6 in spaces 
with positive curvature, e.g., S 2 in Ref. 0. In d = 3, 
these disclinations analogously correspond to bond spin- 
dles. Uncurving the S 3 space that embeds the {3,3,5} 
thus results in irreducible bond spindles with qi.3 > 5. 

Periodic arrangements of 51,3 bond spindles form the 
complex Frank-Kasper crystal structures^, which are 
the densest packings of certain binary hard sphere mix- 
tures^—. Some of these structures even come close 
to approximating the ideal packing bound for Zd (see 
Section IIII A|) . In weakly frustrated fluid configurations 
in d=3, rules for disclination propagation suggest that 
topology constrains the crossing of disclination lines. 
This effect could contribute to the dynamical slowdown, 
because breaking a constraint is an activated event^. In 
this scenario, increasing the fluid density brings the sys- 
tem closer to a perfect tiling of simplices, giving a grow- 
ing importance to these topological constraints, which 
further slows the fluid dynamics. In other words, this 
framework suggests that there should exist a causality 
between the dynamical slowdown and a growing static, 
structural correlation length associated with the bond 
spindle order, which would be an explanation for fragility 
in glass-forming fluids. Note that small frustration is also 
crucial for the avoided critical scenario of the frustration- 
limited domain theory 3 -. 

The small frustration assumption in this argument is 
crucial. Higher-order spindle defects are not necessar- 
ily constrained by the same topological rules, and these 
constraints may go away altogether. Yet regular spaces 
with d > 3 provide no higher-dimensional equivalent to 
the {3, 3, 5}. As detailed in Sectionfl] regular packings of 
simplices in d > 3, are only found on S with a curvature 
that is more than six times that for embedding {3, 3, 5} 
(Table Because even in d—3 the concentration 
of defects is so high that each particle has a multiple de- 
fective bond spindles^, simplex ordering in d > 3 fluids 
should therefore be quite limited. Yet it could be that it 
is not the concentration of defects that matters, but their 
structural correlation, so it remains possible for d = 3 to 
be a limit case, where this geometrical frustration variant 
remains in effect. 



B. Percolating spindle order 

An alternative explanation for the dynamical slow- 
down in d = 3 invokes the percolation of icosahedral 
clusters 43 . These icosahedra are centered on particles 
for which all bond spindles are ideal, i.e., qi^ — 5, which 
result in assemblies of nearly perfect simplices. Once 
these dense and presumably stable structures percolate, 
they have been proposed to affect the dynamical evolu- 
tion of the rest of the system. Because the fraction of 
5i,3 = 5 grows with packing fraction, the onset of icosa- 
hedral percolation should then correspond to the onset 
of sluggish dynamics. And because the percolating struc- 
ture is be expected to become more robust with density, 



this process could also be an explanation for the dynam- 
ical fragility of certain glass-formers. 

Generalizing this approach to higher dimensions en- 
counters two main difficulties. First, there are no equiva- 
lents to the closed-shell assembly of simplices into a gen- 
eralized icosahedra in higher d. One could argue that the 
fraction of qd-2,d = 5 or 4 should be considered instead, 
but then one runs into a second difficulty. Because the 
onset of percolation rapidly decreases with dimensional, 
systems of all densities, even the ideal gas, would be have 
a percolating order. It thus seems that this scenario does 
not apply to higher d, although d—S could still be a limit 
case. 



C. Medium- Range Crystalline Order 

Another variant of the schemes above has been ad- 
vanced by Tanaka and coworker who propose that 
a growing presence of medium-range crystalline order 
(MRCO) underlies the structural origin of the dynam- 
ical slowdown. MRCO would straightforwardly explain 
the dynamical slowdown on the basis that crystal-like 
arrangements are intrinsically dynamically more stable. 
The more of these structural features there are, the slower 
the overall dynamics. Some form of frustration, however, 
must be invoked to explain why MRCO does not extend 
to the whole space, although the fluid is in a domain 
where the crystal is the most stable phase. 

Observations in d — 2 (Ref. |2(J ). and in moderately 
polydisperse systems in d = 3 (Ref. Hi)), have found 
that the growth of MRCO can be correlated with the 
dynamical slowdown. From a dimensional point of view, 
however, this proposal is indefensible. Studies in d > 3 
have found that the fluid order grows increasingly differ- 
ent from the crystal order, and that actually no traces 
of crystal can be found even in monodisperse systems 
where the crystal is unique and singularly dense 3 ^. We 
confirm this result below using a different approach, but 
once again, this description cannot apply to higher d, 
although rf=3 may be a limit case. 

D. Generalized Preferred Cluster 

Coslovich et al. have proposed to generalize the 
above description by systematically investigating the lo- 
cal environment of the particles in a given glass-forming 
fluid to determine the most frequent structural arrange- 
ment o 16 i 17 . This analysis recognizes that the locally pre- 
ferred structure is not necessarily icosahedral or poly- 
tetrahedral, but it also correlates the extension of the 
preferred arrangement to the slowdown of relaxation. 
A related procedure has been put forward by Royall 
et al., with alternative characterizations of the local 
environment of the particles^. It should however be 
stressed that such work focuses on the frequency of a 
given preferred local arrangement, from which one can 
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extract a length scale by dimensional analysis, rather 
than on the spatial correlations associated with the pre- 
ferred local structure. By contrast, a measure of cor- 
relations has been attempted through the use of bond- 
orientational local order parameters, with varying degree 
of succesa^iiiSLSl. 

From a spatial dimension point of view, the notion 
of a preferred cluster is hard to defend. The surface of 
a higher-dimensional hard sphere is so large that none 
of the possible configurations is likely to be significantly 
more represented than the others. But even if a class of 
them were to be more common, this order should be de- 
tected by a Delaunay simplicial decomposition or equiv- 
alently by its dual, the Voronoi decomposition. 



III. GEOMETRICAL MEASURES OF FLUID 
ORDER IN ANY EUCLIDEAN DIMENSION 

In order to determine the extent to which the fluid 
order in d=3 may be singular we consider the develop- 
ment of the generalized spindle order from the Delaunay 
simplicial tessellation in high-dimensional systems. The 
focus is on Euclidean space, as a way to put <i=3 systems 
in perspective. Although the decorrelation property sug- 
gests that all structural correlations in higher dimensions 
can be factorized down to the pair correlation function 
g(r)^, going back and forth between g(r) and the tessel- 
lation is non-trivial. Even the ideal gas geometry (Pois- 
son limit) is not generally solvable. The tessellation also 
offers a differently intuitive description of the packing 
geometry than g[r) alone. 

In order to consistently describe the tessellation or- 
der, we generalize the d = 3 relation established by Cox- 



When d > 3, the analogous formula depends in part on 
the Euler characteristic which is always zero in odd 
dimension and which we take to be zero in even dimension 
as our systems evolve under periodic boundary conditions 
(the periodically fiat torus) . Similar formulae can then be 
derived in all dimensions following the general expression 
we establish in Theorem 151 (Appendix [Bjl . An alternative 
and more numerically stable expression 

Z d = d^. (4) 

1l,d 

is also available in all dimensions (Appendix |Bj . Using 
these general relations, we can obtain two analytically 
tractable limits for disordered packing: the Poisson, i.e., 
ideal gas, limit, and the ideal dense packing limit (Ta- 
ble IIII[) . We consider these two cases before numerically 
evaluating the behavior of hard-sphere fluids at interme- 
diate densities. 




FIG. 1. Schematic of the integration setup for computing 

-ideal 
1i,d 



A. Ideal Limit 

We first consider the densest possible configuration, 
i.e., the ideal packing, following Coxeter's suggestion for 
the statistical honeycomb in Ref. |46|, Sect 22.5. In d=3, 
one imagines the densest configuration obtained by the 
idealized version suggested by Eq. One can wrap 

^ldeai _ __^r_^ regular polyhedra around an edge, and 

thus Zf cai ~ 13.4. The calculation of q[ d % al is straight- 
forward, and provides the ground for generalization: the 
angle between two adjacent faces in a regular tetrahedron 
is arccos(i). This angle can be seen as the central angle 
of a triangle formed with the center of mass of an edge 
of the tetrahedron and the two vertices not on that edge. 
One then divides the total circumference of a circle of 
unit radius (the minimum distance between spheres) by 
the length of the arc cut by this triangle. 

For d > 3, consider a regular d-simplex ad with unit 
side-length and vertices po, . . . , Pd- Pick a z-simplex ai 
in it, say the one formed by the vertices po, . . . , p%. The 
number g^j^ 1 is then obtained by dividing the volume 
of the (d — i — l)-dimensional sphere S centered at the 
centre of mass p of ai and passing through the remain- 
ing vertices by the volume of the section of the sphere 
cut by the non-regular (d— z)-simplex formed by the ver- 
tices p, p»+i, . . . ,Pd- This idea is illustrated in Fig Q] 

The radius r$ of the sphere S satisfies r$ = J (Ap- 
pendix [C]) . It is convenient to rescale the problem and 
obtain gj d ^ al as the quotient of the volume uid-i-i of a 
(d — i — l)-sphere of unit radius centered at by the vol- 
ume of the section of the sphere cut by the non-regular 
(d— i)-simplex formed by the vertices 0, r l5 . . . , r £ ;_ i with 
|rj — ri\ = — whenever j ^ I. Let vi, . . . , Vd-i be unit 
vectors perpendicular to the {d—i — l)-faces of this sim- 
plex touching the origin and facing in the simplex. We 
prove in Appendix \C\ that v, ■ Vj = — 4 . 

Given n vectors vi, . . . , v„ € K™ such that 

[1, U«=Ji 
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we define the region There is a recursive formula for computing V(k,n,p), 

as demonstrated in Appendix [U] 

S(k, n,p) = {x € E™ | ||x|| = 1, and x • v» > p for all i}, 



and let V(k,n 7 p) = Vo\(S(k 7 n,pj). We have established 
that 

-ideal _ ^d-i-1 

9i,c2 — f rj 1 , • n %- W 



Proposition 1 Lei 



f(k,p) = k- 2p 2 k - 2p^l - pVl - k 2 , and 



(n - l)pk + J {I - k) ((n - l)fc + 1) (l - (n - l)p 2 + (n - 2)k) 
* ( *' S (n-2)A= + l ' 



arccos(— k), if n — 2 and p = 0, 

_ Jn/(M,2,0), ifn = 2 andp^O, 

1-3 -,y f k p — kz 



) 2 K -, n — 1, -^^^=^^=^= I dz, otherwise. 



Using this expression, we immediately obtain that 

-ideal _ 27T 
1d-2,d — 7TT- 

arccos(^) 



(7) 



One can also compute the packing occupied by that 
fictitious construction, and obtain, as was done in Ref. |47| 
for d = 3, Roger's bound^. Let 2<7d be a regular simplex 
in M. d of side length 2 having one vertex at the origin in 
R d . Let B d be the ball of radius 1 centered at the origin 
and S' d_1 be its boundary, the sphere in R d . The packing 
fraction (or density) is then obtained by 

_ (d + l)Vol(2cr d n B d ) 
Vol(2<7 d ) ' 

In other words, (fd is the ratio of the volume of the part 
of the simplex covered by the unit spheres centered at its 
(d + 1) vertices to the volume of the whole simplex. 
The packing fraction is thus 

(d + 1)- 



, ddcal 



' Vol(S d - 



2<> 



(d-iy.vd- 



i - , i 

-V(--,d,0), 

which is computed for d=3-8 in Table HU Note that in 
Ref. [28], p. 15, one has that same bound ^ ound on the 

n d/2 ^bound 

center density S, with <pd = — ttt; — • 



B. Poisson Limit 

We now examine the sparsest configuration possible, 
i.e., the ideal gas, which we model by a Poisson process 



TABLE II. Ideal packing fraction ipd in various dimensions 



ideal 
V?3 


— \/2(3 arccos(^) — 7r) 


~ 0.7796 


ideal 
</>6 


~ 0.423(5) 


ideal 
¥i 


~ 0.647(1) 




iff™ 1 


~ 0.33(2) 


ideal 
¥5 


~ 0.5256(1) 




ideal 
¥8 


~ 0.26(1) 



whose geometrical properties self-average in the limit of 
large N. Some of the statistics in which we are inter- 
ested can be found in the encyclopedic Ref. H^. For our 
purpose, we use the results from Ref. [5^ and hence their 
notation. Note, however, that we refer to a "tessellation" 
in lieu of their "mosaic." 

Given X a random Poisson process in R d with intensity 
7, one can find its Voronoi decomposition X and Delau- 
nay decomposition Y. Both are also random processes 
with their own intensities, and their various skeletons all 
have their own intensities. Given a tessellation Z, its j- 
skeleton is the set of all its faces of dimension j. Let 
P^fi be the intensity of the random process and 7^ 
be the intensity of the random process X^>. Because 
the points in X are vertices of the Delaunay decompo- 
sition, = 7, and because the Voronoi and Delaunay 
decompositions are dual to each other, each vertex of Y 
corresponds to a d-dimensional face of X, i.e., 7" = 7. 
More generally, ^ = j( d ~^, from Ref. Hi Thm 10.2.8. 

We can now extend the notion of coordination number 
we defined above. For F S Z^\ let 



N jk (F) = #{S € zW l snz^0} 



and let 



n jk (Z)=E(N jk ) 
be the average value, so rijd — cjjd- 
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TABLE III. Results for Zd- Formulae A and B are formally equivalent, but the former is much more numerically unstable than 
the latter as dimension increases. 
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See appendix for details 





For any stationary random tessellation (the Poisson- 
Voronoi and Poisson-Delaunay tessellations are both sta- 
tionary) of face intensities A*- fc ' , we have (Ref. [50, Thm 
10.1.2) 



\ {j) n jk 



A (fc) 



Tikj. 



(8) 



With certainty, the Poisson-Delaunay tessellation is sim- 
plicial. In a simplicial decomposition, each A:-face is a k- 
simplex hence contains Q+J) j-faces when j < k. Hence 

nkj(Y) — (j^]} when j < k. It turns out that a dual 
phenomenon occur for the Poisson- Voronoi tessellation: 
every fc-face is contained in precisely d — k + 1 top dimen- 
sion cell. Hence nkd(X) = d—k+1 and rijk(X) 
for < j < k < d. We thus obtain 



(d-j+l\ 

fc-j 



-Poisson 
a k,d 



n kd{Y) = j(k) n dk(Y) 



n(d-k),o(X) fd 
no,(d-k)( x ) \k ■ 



v (o) 



(d+l 
\k+l 



( d+1 ) 

\d-k) 



jrn(d-k),o( X ) 



n(d-k),a(X), 



and 



(d-i) 



.(d) 



fl(i) 

= n 01 (Y) = n 10 (Y)^=2 1 

n d , d -i(X) 

= 2 7. HF\ = n d,d~i{X). 

n d -i,d(X) 



To hnd the various rijo(X), we use the relation given 



by Eq. 



the relations 



£(-i)' 



rijk 



k=0 



£(-1)*-%* = 1, 

fc=i 

Et (,) = o, 



Ref.|5fl, Eqn (10.14), 
Ref. [SO, Eqn (10.17), 
Ref.0 Eqn (10.21), 



valid for any stationary random tessellation (including 
the Poisson-Delaunay), the relations 



rijk 



'd - j + 1 
k-j 

fe-i 

(l-(-l) fe )7 W =£(-l 
3=0 



Ref. 



50] p.448, 



from Ref. Ref. [50J, Thm 10.1.5, valid for the Poisson- 
Voronoi tessellation X, and the explicit expression 



v(°) - 



2 d+1 7rl r(^±i)r(i 



d(d + iy. r(f )r(^±i) d r(i) ' 

for the intensity of the vertices of the Poisson- Voronoi 
tessellation (see Ref. |50L Thm 10.2.4). These equations 
contain enough information to fully determine the in- 
tensities 7W, . . . , 7^ and, in turn, all the njk{X) when 
d = 3 (see Ref. [Hfl, Thm 10.2.5 and Table |IV} and when 
d = 4. We can indeed compute 

Theorem 2 Lei X be a Poisson-Voronoi tessellation of 
intensity 7 m R 4 . TTien 
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(0) 
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-7,7 



(i) _ 
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(2) _ 



590 



-7,7 



(3) _ 
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"7; 
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TABLE IV. Limit coordination numbers. 
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177 fy^ 590 



Although this theorem is not explicitly written in 
Ref.— , it was clearly accessible to its authors. We in- 
clude it here for the reader's benefit. 

When d = 5 and d = 6, the same method allows 
one to find the densities 7^ and 7' 1 ). The densities 
j( 2 >, . . . j^ -1 ) are unknown, but are all dependent on 
7^ 2 '. Given that g^-l 8 ™ depends in the missing inten- 
sity, we can use the simulation results to find the missing 
factor (see Appendix [D]), and in turn approximate the 
values of the other quantities. A similar situation occurs 
in d = 7 and d = 8, but two intensity parameters are 
then needed. See Table [TV] for results and Appendix [D] 
for details. 



C. Simulation Approach 

Equilibrated hard sphere fluid configurations with at 
least N = 8000 are obtained in tf=3-8, under peri- 
odic boundary conditions, using a modified version of 
the event-driven molecular dynamics code described in 
Refs. [H and [5l|. Time is expressed in units of y/ flma 2 
for particles of unit mass m at fixed unit inverse tem- 
perature (3. We also consider finite-pressure, hard-sphere 
crystals of D3, D4, D§, and Eg lattice symmetry, which 
form the densest known packings in d=3, 4, 5, and 8, 
respectively. Monodisperse spheres are used, except for 
d=3 fluids, where 50%:50% binary mixtures of diameter 
ratio 7:5 and 6:5 (the larger sphere diameter sets the unit 
length) are used to prevent crystallization, while system- 
atically approaching the monodisperse limit. The densest 
packings of these systems is indeed thought to be frac- 
tionated^,, and the complex phase behavior of similar 
mixtures at finite pressure suggests that the drive to crys- 
tallize is very small^. The properties of these mixtures 
have also been extensively characterized, e.g., Refs. [Tol 
and [Til and Refs. and [H, respectively. 

The diffusivity D is obtained by measuring the 
long-time behavior of the mean-squared displacement 
lim t ^ 00 ((Ar) 2 ) = ^EJr^-r^O)] 2 = 2dDt, while the 
pressure P is mechanically extracted from the collision 
statistics. The QuickHull package 5 - provides a Delaunay 
simplicial tessellation, which is analyzed to extract q~i t d 
and Zd- In d < 6, a tessellation can be directly obtained 
for full periodic configurations, but in d > 6, the tessel- 
lation can only be obtained particle by particle, in order 
to minimize memory usage. In the second approach, the 
values of q^d are obtained by averaging the results over 
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the different particles. In both cases, the Euler charac- 
teristic is explicitly computed to check the consistency 
of the approach. Due to computational limitations, the 
results in d=8 are only averaged over a fraction of the 
particles. 



D. Simulation Results 

As expected, the values of q^d and Zd are found to 
be intermediate between the Poisson and the ideal pack- 
ing limits for all fluid configurations (Fig. ^ . The fluid 
structure evolves smoothly from the Poisson to the ideal 
packing limit, but without ever actually reaching the lat- 
ter. In e?=3, the 6:5 mixture approaches it more than the 
7:5, in agreement with the notion that the former is less 
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FIG. 2. Density evolution of the average bond spindle cjd-2,d 
and number of nearest neighbors Zd with density for fluid 
(squares, dashed line) and crystal phases (circle, solid line). 
The fluid behavior continuously evolves from the Poisson limit 
and reasonably extends to the ideal packing limit (rightmost 
point). The crystal phase, however, clearly sits on a different 
branch and even E& does not converge to ideal limit. (Inset) 
Results for the 7:5 (long-dashed line) and 6:5 (short-dashed 
line) mixtures in d=3, along with the crystal results for a 
monodisperse spheres (solid line). 



geometrically frustrated than the latter. It also allows 
one to reasonably extrapolate that the a non-crystallizing 
monodisperse system would have a qualitatively similar 
behavior that what is found for the mixtures. The fluid 
order thus systematically tends toward an ideal simplex, 
or generalized polytetrahedral, packing. 

As d increases, however, the ideal limit becomes ever 
more distant from the dynamically sluggish regime^. 
Not only the mean, but the entire bond spindle distri- 
bution steadily converges to the Poisson limit (Fig. [3]). 
In that sense, the highest-density equilibrated fluid that 
is dynamically accessible (keeping diffusivity constant) 
gets structurally ever closer to the Poisson limit, rather 
than to an ideal simplex packing at d increases. This 
interpretation is consistent with higher dimensions cor- 
responding to frustrated spaces, because a relatively high 
curvature would be needed to obtain a tiling of simplices. 
Although not a measure of structural correlation directly, 
it suggests that extracting such a length from geometry 
would necessarily shrink with d. Because the static corre- 




FIG. 3. Density evolution of the bond spindle distribution for 
the Poisson limit (thick line), the highest equilibrated fluid 
density attained (solid line - see Fig. |4|, and a couple of in- 
termediate densities. The narrower distribution of the 6:5 
mixture compared to that of the 7:5 mixture at similarly slug- 
gish conditions indicates that the former is closer to the ideal 
packing order the latter. With increasing dimension, however, 
the distribution converges to that of the Poisson limit. 




0P/p 

FIG. 4. The diffusivity decays with increasing reduced pres- 
sure fiP/p in a similar fashion for all dimensions. The be- 
havior of d=3 6:5 (short-dashed line) and 7:5 mixtures (long- 
dashed line) are also similar, although the former initiates 
its slowdown at lower pressures than the latter. (Inset) The 
equation of state for the different fluids. On this scale the 
results for the two d=3 mixtures are indistinguishable. 
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lation length necessary for generating an equivalent slow- 
ing down of the dynamics shrinks with dimension, this 
inference is not in contradiction with static order under- 
lying the process. We note, however, that the order in 
d = 3 systems is not singularly closer to the ideal pack- 
ing limit, even though the ideal simplex packing is much 
closer to Euclidean space than for d > 3. The dynamical 
slowdown in Figure U is also remarkably smoothly evolv- 
ing with d. This continuum of behavior suggests that any 
explanation that provides a singular role to the geometry 
of R 3 cannot reasonably hold. 

Another question of interest concerns the relationship 
between the fluid and crystal branches. For perfect crys- 
tal packings, the crystal's Delaunay tessellation is not 
generally simplicial. Other regular polytopes can be 
identified. Yet when adding thermal noise, and even 
in the limit of infinitely small noise, simplices are re- 
covered with probability one. As a result, in d=3 the 
face-centered-cubic crystal has ^=14, not 12, when ap- 
proaching crystal close packing from lower densities^. 
In higher dimensional crystals, a similar behavior is ob- 
served, but no prediction for the limit behavior of Zd is 
known. In general, the crystal branch is clearly differ- 
ent from the fluid branch, confirming that no crystallite 
form in these systems, in contrast to what happens in 
d=3 monodisperse hard spheres where crystallization is 
facil o 37 i 38 . Even in d=8, where the Eg crystal contains 
an inordinately high fraction of simplices, the fluid order 
follows a branch that is clearly distinct from the crys- 
tal branch. The bond spindle behavior is even located 
outside of the fluid bounds. The dense d=8 fluid thus 
bears no trace of the high-simplicial nature of the crys- 
talline phase. Actually, even the crystal structure is but 
partially captured by the ideal packing. Although the 
behavior of Z is similar to that of the other lattices, the 
bond spindle behavior is not even bounded by the ideal 
packing. We stress that no MRCO is found in high-rf 
fluids, confirming the conclusions of earlier fluid order 
studies^. 



IV. "ORDER AGNOSTIC" LENGTHS 

The results of the above analysis indicate that the 
simple geometrical considerations brought forward to ex- 
plain the dynamical slowdown are unlikely to provide any 
significantly growing static correlation length in hard- 
sphere glass-formers under compression. Other mea- 
sures of such a correlation length, e.g. through bond- 
orientational order, confirm this conclusion (see Ref. [27| 
for the case d = 3) . Yet this analysis does not exhaust the 
possibilities for more subtle forms of structural correla- 
tions. Other proposals have been made based on nontriv- 
ial "order agnostic" lengths^, such as the patch repeti- 
tion length associated with the notion of configurational 
cntrop y 55 i 56 , length scales extracted from finite-siz o 57 i 58 
and information theoretic analysis^, as well as point-to- 
set correlation length o 23 ' 58 ' 60 " — . 



The point-to-set correlation lengths, in particular, can 
be studied by considering the distance over which bound- 
ary conditions imposed by pinning particles in a fluid con- 
figuration affect the equilibrium structure of the remain- 
ing (unpinned) particles. The original proposal, moti- 
vated by the random first-order transition theory^, con- 
sidered a cavity whose exterior is a frozen fluid configu- 
ration^. The crossover length characterizing the range 
over which the pinned boundary determines the con- 
figuration inside such a cavity also enters in a formula 
bounding from above the relaxation time of the (bulk) 
fluid 23 . Yet other geometries of the set of pinned par- 
ticles also allow one to extract point-to-set correlation 
lengths—. These lengths need not coincide nor evolve in 
exactly the same way as the system becomes sluggish, 
but they nonetheless characterize the same underlying 
physics. Evidence for the growth of such lengths as the 
relaxation time increases has been found in several model 
glass- formers^— . 

These point-to-set correlation lengths are more general 
than structural lengths based on a specific local order de- 
scription and are thus expected to provide upper bounds 
for the latter. Although they may provide less detailed 
geometrical information, and may be numerically harder 
to investigate over a broad dynamical range, it seems 
clear that a correlation between structure and dynamics, 
if present, should at least emerge at their level. For this 
reason, we complete our study of structural lengths re- 
lated to ideal simplex order in hard-sphere glass-forming 
fluids with that of the "order-agnostic" point-to-set cor- 
relation length obtained from a random pinning of a frac- 
tion of the particles in equilibrium configurations. 



A. Overlap Function 

We consider a situation in which a fraction c of the 
particles of an equilibrium fluid configuration is pinned at 
random. Information about a possible static length scale 
growing with the dynamical slowdown as one compresses 
the fluid is then obtained from the long-time limit of 
the overlap between the original configuration and the 
configuration equilibrated in the presence of the pinned 
particles. If the reference and the final configurations 
are quite similar, the average pinning spacing is shorter 
than the static correlation length, and the opposite if 
the two configurations are dissimilar. Previous studies 
have chosen an ad hoc boxing of spac o 27 ' 61 " — , but in 
systems where the particle density is changed, it is more 
robust to measure the configuration overlap by means of 
a microscopic overlap function 

w mn (t) = e(a-\r n (t)-r m (0)\), (9) 

where a = 0.3a is chosen sufficiently small to enforce 
single overlap occupancy for hard spheres. For a concen- 
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tration c of pinned particles we therefore have 



Qdt) 




(l-c) 2 N 



(10) 



where the brackets denote an average over equilibrium 
configurations, the over line represents an average over 
the different ways to pin a fraction c of the particles of 
a given equilibrium configuration, and the sum is over 
all unpinned particles, with B denoting the set of pinned 
particles. 

When t = 0, pinning is irrelevant and the overlap 
Qc(0) = 1. When there are no pinned particles, i.e., 
for c = 0, 



h(t) 



1 

N 



dr J dr'e(a-|r-r'|)(p(r,0)p(r',t)), 

(11) 

where p(r,t) — Yl n =i ~ r «W) i s tne density of par- 
ticles in a configuration {r^i)} at time t. The overlap 
Qo(t) could be expressed in terms of van Hove functions 
whose weight is provided by the microscopic overlap func- 
tion, but this approach would not be not particularly il- 
luminating. Instead, at t — 0, we use the equilibrium 



,66 



property 



p5(r-r') + p 2 [l + p 2 h{\r-r'\)}, (12) 



where h(r) = g(r) — 1 is the total correlation function, 
together with the fact that for hard spheres h{r < a) = 
— 1 (and thus also when particles are within a distance 
a < a), to recover that Qo(0) = 1. When t —> oo, one 
finds 



Qo(oo) = ^ J dr J dr'e(a- 



-v'\)p 2 = P V d {a). (13) 



The relevant quantity one wishes to examine is thus 
Q c (oo) — Qo(ao). By locating the rapid growth of the 
overlap with c~ 1 / d , we can extract the crossover length 
between small and large overlap, which corresponds to 
a static point-to-set correlation length within the bulk 
system. 




FIG. 5. Growth of the overlap with increasing pinning con- 
centration in d = 3 for two different hard-sphere mixtures. As 
the system density increases, the crossover from low overlap 
to high overlap takes place at an ever smaller concentration 
of pinned particles. The average spacing between defects at 
a 0.4 overlap defines £ p . (Inset) Time evolution of the over- 
lap and of the self component Q s c . The long time value of the 
overlap is attained when the self-part has completely decayed. 
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FIG. 6. The static (point-to-set) length £ p for d=3 6:5 (short- 
dashed line) and 7:5 (long-dashed line), extracted from the 
data in Fig. [5] both saturate the bound from Eq. (|15[) (thinner 
lines). In higher dimensions, the bound becomes continuously 
weaker. The £dyn results for the 7:5 mixture (dot-dashed line) 
are taken from Ref. [TlL 



B. Simulation approach and results 

We have calculated the point-to-set static length from 
the pinning analysis in d=3 only. In higher dimension, 
the computation becomes prohibitively time consuming. 
Also, smaller systems with = 1236 are used and the 
absence of finite-size effects are verified for select systems 
with N = 9888. The length £ p is defined as the value 
of the mean distance c~ x / d between pinned particles for 
which the overlap falls below 0.4 (Fig. [5]). The extracted 
length is not very sensitive to this choice, provided it is 
intermediate between low and high overlap. 

Figure [6] shows that the static length £ p increases only 
very modestly over the dynamically accessible density 



range, while the diffusivity D and the structural relax- 
ation time r a change by about 4 orders of magnitude 
(Fig. @|. Interestingly, over the same density range, the 
"dynamical length" £dyn characterizing the extent of dy- 
namical heterogeneity grows markedly (the results shown 
for the 7:5 mixture are from Ref. IllT l. Whereas the 
change in the static length is measured in fractions of 
their low-density value, £dyn grows by a factor of 4.5 when 
reaching f3p/ p — 22, and seems to keep on growing^. 

We can verify that the static length £ p extracted sat- 
isfies the bound derived by Montanari and Semcrjian 23 
that relates the structural relaxation time scale and the 
(cavity) point-to-set length £ps- One may indeed rea- 
sonably expect that £ p ~ £pg, as long as one is far from 
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any putative random first-order transition 67 , which is the 
case here. (When approaching a random first-order tran- 
sition, one should have that £ p ~ Cps-) The bound states 
that 



< 



TO 



cxp (Bg s ), 



(14) 



where tq is a constant setting the microscopic time scale. 
The coefficient B depends on pressure (and tempera- 
ture, more generally) and is such that when £ps ~ a 
the right-hand side describes the "noncooperative dy- 
namics" of the model 6 ^. Using an Arrhenius-like ar- 
gument for activation volumes 1 —, we note that, all else 
being equal, higher pressures trivially rescale the free- 
energy landscape and thereby slow the dynamics. For a 
hard-sphere fluid, one then expects B oc /3P. The upper- 
bound of r a diverges with the pressure even in the ab- 
sence of any growing £pg , as when approaching T = for 
an Arrhenius temperature dependence. In the low and 
moderate density fluid, the relaxation time indeed fol- 
lows T a (P) ~ riow(-P) = T oexp[X/3P] with K a density- 
independent constant. One then finds that 



frs(P) > ( log[r a (P)/r ] 

(P)/M 



l/d 



(15) 



where £ps,o is the low-density limit of £ps- Equation (|15[) 
thus provides a lower bound for the growth of a static 
length imposed by the dynamical slowdown. As seen in 
Fig. [6l the bound given by Eq. (|15j) increases only slowly 
in the dynamically accessible domain and is already sat- 
isfied by £ P . 

In dimension higher than 3, we are unable to directly 
compute the point-to-set length £ p , but we can nonethe- 
less use the bound in Eq. (flB"]) as a proxy for the growth 
of a static length. The results in Fig. [6] give strong evi- 
dence that the growth of any static length is very limited 
for hard-sphere glass-formers in the accessible dynamical 
range where the relaxation time increases by 3 to 4 orders 
of magnitude. 



C. Linear response for the overlap at small 
concentration of pinned particles 

We conclude this section on order-agnostic static 
lengths by proving that one cannot gain any knowl- 
edge about point-to-set correlations from a mere inves- 
tigation of the small-concentration pinning regime. In 
this limit, Q c (oo) — Qq(oo) indeed only conveys infor- 
mation about the equilibrium fluid pair correlation func- 
tion. Higher-order terms in the expansion in powers of 
c involve higher-order equilibrium correlation functions 
of the fluid, the order p of the expansion being typically 
expressible in terms of n-body correlation functions with 
n < p. 

The overlap between the original equilibrium config- 
uration and that equilibrated in the presence of pinned 



particles after due average can be expressed in the long- 
time limit as 



dr f dr'e(a-|r-r'|)(p(r) (p(r')) e ) 



Q c (oo) = 
1 

(1 -c) 2 N 

(16) 

where p(r) = J2 n ^ B S(r — r„) is the local density of un- 
pinned particles and ( ) B represents a conditional average 
over configurations of particles at equilibrium with a set 
B of pinned particles. Specifically, 



(p(r')) B = / dr 



e -/3V({f N }) 

Z N ({r m },B) 



l[5(v k -r k )p(r'), 

(17) 



keB 



where 

4({rJ, B) = / dr N e -^ v « fN » [J 6(v k - f k ), (18) 

J keB 

V({r N }) denote the interaction potential, which corre- 
sponds to pairwise hard-core terms here, and p(r') = 
Sm£B l K r ' ~ *m)- One should keep in mind that (p(r')) B 
still depends on the reference configuration {r m } and on 
the pinned set £>, the integration being on the configura- 
tions denoted {f^} = {r m }. 

The conditional probability that is used in Eq. p7|) 
satisfies the requirement that when integrated with re- 
spect to the equilibrium reference configuration it gives 
back the unconstrained equilibrium distribution 6 -^—. As 
an illustration, one obtains 



({p(r')) B ) = / dr 



-/3V({r N }) 

Zn 



Y[S(r k -r k )p{r')= / dr N 



dr 



keB 



Zn 



-/3V({f N }) 

Z N {{r m },B) 

'}) 

— p(r')x 



dr 



N 



-/3V({r N }) 



Z N ({i m },B) 



-r fe = 



keB 



dr 



a -/3V({r N }) 



■ P(V) = (p(r')) , 



where use has been made of the delta functions to change 
the arguments of ZN({r m }, B) and of the definition of the 
constrained partition function in Eq. (|18p. 

To organize the expansion in powers of the concentra- 
tion c, rather than considering realizations with exactly 
cN pinned particles, it is easier to consider the situa- 
tion in which cN particles are pinned on average, as in 
a grand-canonical treatment- 7 - -^- . Let r be the set of TV 
occupation variables {t^ = 0,1} characterizing the spe- 
cific pinning of particles in a given configuration. The 
corresponding probability distribution is then simply 



P n (t) = (1 - c) 



(20) 
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We then rewrite the two-point density appearing in the 
expression of the overlap Q c (oo) as 



p(r)(p(r')) g =tr T P N (r) / dr 



-/3V({r N }) 



-p(t;t)x 



7 Zjv({r m },B(r)) .A* _, J 



feeB(r) 



(21) 

where we have now explicitly indicated the dependence 
on the pinned set r. 

The desired expansion in concentration c can be gen- 
erated by first formally expanding the average over the 
pinned set in increasing number of pinned particles. For 
a generic function J"(r), this reads 

f 1 N N 

J"(r)=trJ ^- C p(1-c) n -p J2 11^,1) 
^ p>0 ki,..,k p =l i=l 

x [] Jfo,0) ^(t) } , 

(22) 

where 5(a, &) is the Kronecker symbol. After expanding 
the remaining factors in (1 — c), one finally arrives at 

r N N 



F(t) =tr 7 



j=i 



k=l 



JV 



Il^ ) + T E (S(r kl ,l)-6(r kl ,0))- 



&k 



fei ,fc 2 = l 



(«5(r fe2 ,l)-«J(r fe3 ,0)) [] <5(t„0) 



(23) 

where only the first terms have been explicitly given. 

Eq. fl23J| should now be applied to (p(r) (p(r')) e ). The 
zcroth-order term corresponds to the case with no pin- 
ning (B is empty and p(r) = p(r)) and simply gives 
(p(r) (p(r'))) = p 2 . The linear term can be expressed 
as 



E E /* 

fc— 1 m.n^k 



S(r-r m ) dr N * 



<J(r fc - r fe )<S(r' - ?„) - JV V) / dr 

/ 



Z]v({r m },fc) 

e -/3V({r N }) 



the fluid in the absence of pinning bb 

/' e -/3V({r N }) 

iV / II dr *^ = p«(r fe )=p, 



N(N-l) H dr, -— =p WQ Tk - Tkl \) 



i^k,k' 

one can rewrite Eq. (f24|) as 



(26) 



iV(iV 



7J2 E J dTl J dr ™ J dr n 5{r -r m )8(r' -v n ) 



x p (2 )(|ri - r 



\) p W(\ ri -r n \)-Np\ 



(27) 

Because a cancellation of the leading terms in N is an- 
ticipated, one has to be cautious about the sub-extensive 
terms in the long-distance limit of the pair density, i.e. 
p( 2 >(r) — ¥ p 2 (l — a/N) when r — > oo. (For an ideal 
gas, p (2) (r) 7V(jV - l)/V 2 so that a = 1.) In con- 
sequence, we define p( 2 '(r) = pS(r) + p 2 [l — a/N + h(r)] 
with h(r -» oo) = + £>(l/iV 2 ) - compare with Eq. ([T2]>. 
Noting that the <5 terms in the pair densities do not con- 
tribute, Eq. (f2"T|) can then be expressed as 

-2p 2 a-pjdvh{r) + p 3 J dr"h(\r" -r\)h{\r" -r'\) 

(28) 

where sub-extensive corrections of 0(1/ N) have been 
dropped. We now determine the coefficient a by requir- 
ing that (p(r) (p(r')) e ) -¥ (1 — c) 2 p 2 when |r - r'| -> oo. 
When expanded in powers of c, this condition implies 
that a — p J drh(r) = 1 (as indeed found for the ideal 
gas). 

After inserting the above results for the zeroth and first 
orders in Eq. (|16[) . subtracting the value of Qo(oo), and 
expanding the term in (1 — c) 2 in the denominator of the 
right-hand side of Eq. (|TB)) . one finally obtains 

Q c (oo) - Qo(oo) = 

cp 2 / / drdr' e(a-r)h(r')h(\r' - r|) +C(c 2 ) . ^ 



-8{r - r m ) 



„-/3V({r N }) 

dr^^— <5(r'-f„ 



(24) 



where because of the fluid's translational invariance 

Z N ({r m }, k)=J d~r N e-^« fN »<5(r k - f k ) = ^ . 

(25) 

By using the definition of the 1- and 2-body densities of 



For an ideal gas (h = 0), one recovers the expected result, 
Q c (oo) - Qo(oo) = 0. 

We have focused here on the linear term in c which ex- 
presses the linear response of the fluid to a perturbation 
associated with pinning a vanishingly small concentra- 
tion of particles. The higher-order terms in powers of c 
can be similarly computed, but lead to increasingly te- 
dious algebra. It is easily realized that the pth order then 
involves up to p-point density correlation functions. 

As is clear from the above expression, the linear co- 
efficient of the dependence of the configuration overlap 
on the concentration of pinned particles does not carry 
any relevant information on point-to-set correlations. It 
only involves the pair density correlations. As a result, 
one can at best hope to extract from it the correspond- 
ing two-point density correlation length, a quantity that 
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is not of much interest in the context of glass-forming 
fluids. 

In AppendixjEj we also relate the configuration overlap 
to the "blocking" (or "disconnected") correlation func- 
tion introduced in the statistical mechanics of fluids in 
a disordered porous medium, and more generally in the 
theory of disordered systems. 



case, the vertex figure {q} contains obviously #Fp = q 
points and #Fp = q edges, corresponding to #-F 2 = 
#F 3 = q hence g ,2 = «?• □ 

Appendix B: relationship between Z and 
coordination numbers 



V. CONCLUSION 

In summary, our results demonstrate that traditional 
descriptions of geometrical frustration are of limited ap- 
plicability to hard-sphere and related systems, and that, 
more generally, structural correlations are fairly limited 
in the dynamical regime accessible in simulations. The 
complex and related question of the relationship between 
the structural and the dynamical correlation lengths will 
be the subject of a future publication. 
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Appendix A: Simplex Coordination and Schlafli 
Notation 



Let fj be the number of aj in a ad, i.e., the number of 
faces of dimension j. 

Lemma 4 We have fjNd — q~jNj, for < j < d and 
Ni = NZ/2. 

Proof: The second relation is an immediate conse- 
quence of the classical hand shake lemma: in a graph, 
twice the number of edges must equal the sum of the de- 
grees of the vertex, so 2N\ = ~}2,z(i) = NZ. Similarly, 
if each ad in the decomposition has a number fj of faces 
aj, we must have fjNd — X)j = QjNj- □ 

Theorem 5 Let x oe the Euler characteristic of the d- 
dimensional manifold in which we compute the coordina- 
tion numbers. Then 



Z,, 



d(d + l)(l- X /N) 



Note: The Euler characteristic is automatically zero 
when d is odd or when the manifold is a flat torus. As 
discussed earlier, the Euler characteristic is always zero 
in this paper 

Proof: We have 



Lemma 3 Let {ai,...,a n } be a regular honeycomb in 
W l or S n . Then, at each vertex, q n -2,n = a n- 

Proof: Let p be a vertex of the honeycomb and F be 
the subset of all the faces in the skeleton of the honey- 
comb that include p but are not p. One has the map- 
ping between F and the skeleton of the vertex figure of 
p, sending each <7j to a <7j_i. Let tt denote this map. It 
is an inclusion preserving bijection. Let F p be the full 
skeleton of the vertex figure at p. Filter F and F p by the 
dimension of the cells by an upper index. One therefore 
has that ir restricts to a bijection F l — > F^ 1 . 

Let P be the honeycomb and P p be the vertex figure 
of P at p. Given s £ F l , we have 

q hd { Sl P) = #{teF d \s(lt} 

= #{i G Fp | s C t} (since p G s anyways) 
= e F^ | ttOO C n(t)} 

= ?i_i,d-i(7r(s),P p ). 

Hence one can prove this claim by induction. One goes 
down to the lowest dimension where qd-2 makes sense, 
that is when d — 2. In this case, qd-2.d is the number 
of polygons around a particular vertex. In the Euclidean 



x = Y,(-l) j Nj 





-JVi-f 




= N — 


NZ 
~2~ ' 


3=2 * 


= N — 


NZ 
~2~ ' 




= N — 


NZ 




Note that fj = 




, in particular fo = d + 1 



fd = 1. Note also that, obviously qd = 1, and qd-\ = 2. 
Solving, we get 



Z = 



2(1 - x/N) 



3=2 \ 2 ) 

which simplifies to the desired formula. 



□ 
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Note that qo doesn't appear in this formula. An alter- 
native formula is 



Z = d 



(Bl) 



Indeed, repeated use of lemma 0] tells us that 



fo fo qiNi f qiZ 



fd+l\ 



N h h 2 



Appendix C: ideal coordination numbers 

In this appendix, we cover two aspects of the compu- 
tation. First, we prove Proposition [TJ Second, we prove 
the claims left unproved in the establishment of Eqn (|5|) . 



1. Proof of Proposition Q] 

It is useful for the proof to introduce one more variable: 
the radius of the sphere. Given n vectors vi, . . . , v„ G M" 
such that 

j k, if i ^ j, 



1, if i 



we define the regions 



T(k, n,p) = {x G R" | x • Vi > p for all i}, 
S(R, k, n,p) = {x G M™ | ||x|| = R, and x G T(k, n,p)}, 



and let 



V(R,k,n,p) = Vo\(S(R,k,n,p)). 



To conform with the notation in the main text, we have 
V(k,n,p) = V(l,k,n,p). 

Lemma 6 We have V(R,k,n,p) = R n ^V(l, k,n, -^). 

Proof: If x G S n ~ 1 (R) satisfies x • Vj > p then G 
satisfies J| • Vj > ^r, and vice- versa. So the 
map x i — y gives a bijection between S(R,k,n,p) and 
S(l,k,n, J^). This map shrinks hypersurface volumes by 

R"- 1 . a 

We can arbitrarily rotate the v^. So let us assume that 



' n — c n 



(0,...,0,1). 



x K, we 



Lemma 7 Suppose v„ — e n . For (y, z) G R r ' 

have 

z > p, and 

(y, z) G T(fc, ra,p) T ,_fc_ _ 1 _p-fc^x 



Proof: Since v n = e„, we have that ■ e n — k for all 
i < Let w, = Z z Ztt"u , so that w, G M n_1 . Since 



Ivj - fcejl 2 = 1 — 2A: 2 



fce„. Moreover, for i ^ j, we have fc = • Vj = (1 — 
k 2 )(\Vi,Wj) + k 2 , hence 



Wj • W.; = 



k- k 2 
T^k 2 



k 



1 



Note that for i < n, x • = (y + ze n ) ■ — /c 2 Wi + 
ke n ) = vl — k 2 y-Wi+kz. Hence the conditions x-v* >p 
is equivalent to 

p-kz 

y.w i >- 7 _,for*<n. 

The condition x • v n p is equivalent to z ^ p. EH 

Lemma 8 Suppose v n = e„. TTien (y,z) G S(k,n,p) if 
and only if 



p < z < h(k,n,p), and 



p — kz s 

VT - fc 2 "' 



Proof: Let's start with the easy part. Let x = (y, z) G 
IR"" 1 x R. If |x| = 1, then |y| = Vl - z 2 . It is also 
clear that we must have z > p, at least from the previous 
lemma. 

Geometrically, we see that for (y, z) to be in that re- 
gion, the condition z > p is not enough. The "tip" of the 
region occurs at a point X where X- Vj = p for 1 < i < n, 
and X • e„ = h, and |X| = 1. 

Let g be the matrix with entries gij = Vj • Vj. Let X 
be the column vector of the coordinates of X in the basis 
Vi. So X = J2 x * v i- Let Y = [p ••• p h]. Then 



Yi = X • v, = ^ XjVj • Vi = ^ ftXj. 



so Y = g±, or X = g^Y. Then 

1 = X ■ X = %%9ij = X T gX = Y T g _1 Y. 

ho 

So h is a solution to the quadratic equation Y T g~ 1 Y = 1. 
There are two solutions to this equation, we need to take 
the solution greater than p. The function h = h(k 7 n,p) 
is that solution. □ 

We have now set up the stage to perform a volume com- 
putation in cylindrical coordinates. One has to be quite 
careful in using cylindrical coordinates on the sphere. 



Lemma 9 We have the relation 
R 



dvo\ sn{R) 



VR 2 



dvolo„-i 



dz 



between the volume element on the sphere of radius R 
centred at the origin o/R™ +1 and the volume element on 
the parallel sphere at height z. 
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Proof: Let s n : W 1 — > 1" be the inverse of the stere- 
ographic projection, and 

q: W 1 - 1 x [-R,R] ->■ R n+1 

(x , t) i y {VP 2 -t 2 s n (x),t). 

The image of q is all of S n (R) less half of a great circle. 

Recall that (s^gs^-^x — (ds n ) x (ds n ) x . Similarly, we 
have that (q*gs^(R)) = (dq) T (dq). But 



(dq)( x ,t) = 
so the metric is 
g = {dqf \dq) = 



VR 2 - t 2 {ds n ) x VR2 _ t2 s n {x) 



(R 2 ~t 2 )s* n g Sn - 1{1) 

(R 2 - i 2 )< 5s „-i (1) 




-t(dSn)x S n {x) 



R 2 -t 2 



R 2 



1 



R 2 ~t 2 



By taking the determinant and a square root, we get 



dvol 



S»(H) 



R(R 2 - t 2 ) — dvol S n-i W dt, 



as desired. Rescaling, we establish Eqn 



□ 



Proposition 10 For n > 2, we have 



V(k,n,p) 



h(k,n,p) 



{l-z 2 )^V(- 



+ k 



n — 1, 



p — kz 



VT^WVT^T 2 



)dz. 



Proof: Let f\ 1 f2, fj,{-,z)he the characteristic functions 
of the set S(k, n,p), the interval [p, h(k, n,p)] and the set 
S(\/l-z 2 , j^,n-l, ^=f) respectively. So /i(y,z) = 
f2{z)fs{y,z) by Lemma [8] We have 



V(k, n,p) 



fi(y,z)dvol S n-i(y,z) 



s n - 
l 



VT^ 2 



1 JS"- 2 
h(k,n,p) V(y/l -Z 2 , j^, 



rfvol s „_ 2(VT ^ ) (y)dz 



1,4 



dz. 



VT^z 2 

Then, using Lemma [SJ we obtain the desired result. □ 

To conclude the proof of Proposition [1] we need the 
following two lemmas. 

Lemma 11 V(k, 2,0) = arccos(— k). 

Lemma 12 V(k,2,p) = V(f(k,p), 2, 0). 

Proof: Figure [7] illustrates the idea of this proof. Let 
ni, 112 be unit vectors perpendicular to vi and V2 re- 
spectively, pointing in the direction of the region we are 
computing. Let fi! to be the vector closest to ni on the 
circle S 1 and in the plane x • vi = p. Let A2 be similarly 
defined. Let Vi and V2 be the vectors such that 



x • Vj > p 



x • v, > 0. 




FIG. 7. Illustration for the proof of Lemma [T2l 



We have v 4 • = 0. Let k* = \J\ -k 2 . 

Let a = angle(vi,V2) and 9 = angle(ni, 112). We have 
a = 7r — 9. Thus 



ni n 2 

and 



cos(#) = cos(7r — a) = — cos(#) = — vi-V2 



IT 7T 

112 • vi = ni • V2 = cos(— — a) = cos(# — — ) = sin(6') = k* . 

We have ni = p*ni + p\~i and n.2 — p*ii2 + pv2- Thus 
ni • n 2 = (p*) 2 ni ■ n 2 + 2p*pvi ■ n 2 +p 2 vi • v 2 



1G 



-(p*) 2 k + 2p*pk* +p 2 k 

Hp 2 - (p*) 2 ) + 2 P *pk* 

k(p 2 - (1 -p 2 )) + 2p*pk* 



= 2p 2 k -k + 2p\J 1 — p 2 \J\ - k 2 
= -f(k,p). 

And thus Vi • V2 = f(k,p). □ 

2. Proof of Eqn. ([5]) 

We know that po, . . . , Pd are equidistant, with \pj — 
P/| = 1 for j < I. Let pj = pj — po. Then for j ^ I, we 
have 1 = |pj - p ; | 2 
hence 



P]-Pi\ 2 = IPjI 2 + IPil 2 — 2pj • Pz 



Pj ■ Pi 



The centre of mass of <Ji is p = p "\'^ Pi and the distance 
square from p to p n +i is 



|p- P0 | 2 = |p-(i + l) P ° 12 



i + V 

(Po - Po) H h (Pi - Po) 12 



t+1 

riw E Pi -pi 



(i + 1) 2 

1 r i(i-l) 

(i + iy K 2 

i 



+ 



2(2+1) 



For j = (i + 1), . . . , d, the triangle PoPPj has a right 
angle at p. Hence we must have Vi = |pj — p| = 

y/i - |p - p| 2 = ^1 - = v ■ 

Let Conv(p i+ i - p, . . . , p d _i - p) and Conv(p l+2 - 
p, . . . , Pd — p) be two (d — i — l)-dimensional faces of 
the (d — z)-dirnensional simplex that is the convex hull of 
p, Pi+i, . . . , Pd- Let vi and v 2 be unit vectors perpen- 
dicular to these faces, facing in the simplex. 

For simplicity of notation, let's relabel the vectors: let 
Pj = Pi+j — p and m = d — i. So now our vertices are 
pi, ... , Pm and 0. Let Sj = 21+2- We have ||pj || = r\ and 
for j I, since 1 = ||pj-p/|| 2 = ||pj || 2 + ||p/|| 2 - 2p r p/, 



Pj Pi 



Using the analog of the cross-product in higher dimen- 
sion, we let 

wi = *pi A • • • A p m -i, 

w 2 = (-^♦pjA-Apm. 



Then Vj = Wj/|jwj|j if the vector pi, . . . , p m form a basis 
compatible with the orientation. We have 



Wi|| 2 = || W 2 || 2 = dct(p, • Pl)2<j,l<r 



det(S jir 2 + (l-*ii)*i) 2 <j,i<m 
(r 2 + (m-2) Sl )(r 2 -s l ) m - 2 , 



while 



w r w 2 = (-1)™ 1 det(p j -p i+ i)i< i j< m _i = -{rf-s, 



\m— 2 



Sj. 



So 



Vl • v 2 



Wi • w 2 
|wi||||w 2 | 

Si 



rf + (m — 2)sj 



1 



Appendix D: Statistics for Poisson Voronoi 
tessellations in high dimension 

1. Dimension 5 

Theorem 13 Let X be a Poisson-Voronoi tessellation 
of intensity 7 in R 5 . Then 

(m 1296000 4 m 3888000 4 

7 V ' = 7T 7, 7 V ; = 7T 7, 

' 676039 '' ' 676039 ' 
(2) 3888000 4 (4) 

■-yV ; = j[ /y _L ry\ I ry 

' 676039 ' ' 
,„ 1296000 4 / 4 \ 

7 (3) = ^4 + 2 (4) _ 2 

' 676039 ' ' ' 



-Poisson 
93,5 — 


19440000tt 4 7 




3888000tt 4 7 + 676039 7 ( 4 ) 


- 6760397 ' 


-Poisson 
92,5 — 


12960000tt 4 7 




648000tt 4 7 + 676039 7 ( 4 ) - 


676039 7 ' 


-Poisson 

9l,5 - 


19440000 tt 4 7 
676039 7W 




-Poisson 
90,5 


7776000 4 
676039 W ' 




^Poisson 
Z 5 — 


2 7 W 

7 





while in all other cases where j < k, we have 

nk ^ = ^m{k-j 

Numerical simulations provide an estimate for the 
value of 53 ° lsson . Using this estimate, we find an esti- 
mate 7 (4) ~ 44.206608727. Substituting this value in the 
quantities above yield estimates for the other coordina- 
tion statistics agreeing with the simulated value with a 
relative error of 0.5%. By opposition, using the known 
value of g 3 d | al to perform the same trick yields quanti- 
ties agreeing with relative errors 16% for 9i* 5 lsson , 13% 
for Z% oisson , and 53% for Zf cal . 
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2. Dimension 6 

Theorem 14 Let X be a Poisson- Voronoi tessellation 
of intensity 7 in R 6 . Then 



7 



(o) 



12964479 



7 



0) 



10000 
5 7 ( 2 ) 90751353 



-. 7 « = 

7, 



90751353 
20000 



-7, 



2 8000 
m „ r ,N 220356143 

7 (4) = 2 7 ( 2 ) 7, 



7 



(2) 



20000 
116560311 
40000 



-Poisson 
94,6 — 


272254059 7 
100007( 2 ) ' 


-Poisson 
93,6 


3630054127 


200007( 2 ) - 90751353 7 ' 


-Poisson 
<?2,6 — 


907513530 7 


400007(2) - 2203561437' 


-Poisson 

Qi,e — 


1089016236 7 


200007(2) - 1165603117' 


-Poisson 

%,6 — 


90751353 
10000 ' 


^Poisson 
Z 6 


200007(2) - 1165603117 


2OOOO7 



while in all other cases where j < k, we have 
n kj (X) 



7 0) (7 - J 
7W \k - j 



Numerical simulations provide an estimate for the 
value of ^4 g lsson . Using this estimate, we find an esti- 
mate 7(2) ~ 6030.5244987. Substituting this value in the 
quantities above yield estimates for the other coordina- 
tion statistics agreeing with the simulated value with a 
relative error of 0.3%. By opposition, using the known 
value of g 4 d 6 al to perform the same trick yields quanti- 
ties agreeing with relative errors 25% for qi™ sson , 20% 
for Zf oisson . 



3. Dimension 7 

Theorem 15 Let X be a Poisson-Voronoi tessellation 
of intensity 7 in R 7 . Then 

m 3442073600000 6 

7 — n 7' 



7 M = 



v( 3 ) 



322476036831 
137 



7T 6 7, 



v(4) 



322476036831 

48189030400000 „ „ m 

7rS + 3^ ( ) 

322476036831 ' 

58515251200000 



7 ( 5 ) = 



-Poisson 
#5,7 



-Poisson 
<?4.7 



107492012277 



7T 6 7 + 3 7 ( 2) + 7 ( 6 ) - 7, 

tt 6 7 + 7(2) + 2 7 ( 6 ) - 2 7 , 



96378060800000 tt 6 7 
322476036831 7(2!' 

192756121600000tt 6 7 
9674281104937(2) - 48189030400000tt 6 7 ' 



-Poisson 
93,7 



-Poisson 
<?2,7 



240945152000000tt^ 



7 



9674281104937(2) + 32247603683l7( 6 ) - (58515251200000tt 6 + 322476036831) 7 ' 

192756121600000tt 6 7 
3(1074920122777(2) + 2149840245547( 6 ) - (6884147200000rr 6 + 214984024554) 7 ) ' 



-Poisson 

yi,7 



-Poisson 
<?0,7 



^Poisson 



96378060800000 tt 6 7 
322476036831 7^' 
27536588800000 6 
322476036831 n ' 

2 7 ( 6 ) 

7 

while in all other cases where j < k, we have 

n ^ x ^^m{k-j 

Numerical simulations provide an estimate for the 
value of q^° 7 mson and g|^ lsson . Using these estimates, 



we find an estimate 7 (2 ) ~ 64873.075327 and 7(2) ~ 
228.84946857. Substituting these values in the quantities 
above yield estimates for the other coordination statistics 
agreeing with the simulated value with a relative error of 
0.3%. 



18 



4. Dimension 8 

Theorem 16 Let X be a Poisson-Voronoi tessellation 
of intensity 7 in M 8 . Then 



v( 7 ) - 



68567707769134446 9 7 ( 2 
45956640625 7 ~ 



v (4) 



7 



(o) _ 



v(!) = 



4155610296921974 

45956640625 
18700246336148883 



7, 



7, 



v(3) - 



v(5) - 



v(6) - 



45956640625 
12466830890765922 
6565234375 
37400492672297766 



-7 + 



7 7 



(2) 



7' 



6565234375 
257647930322443638 
45956640625 



4 1 
17 7 (2) 



57 



(4) 



2 ' 
2 7 (4) , 



-Poisson 
^6,8 



-Poisson 
?5,8 



-Poisson 
94,8 



149601970689191064 7 



45956640625 



7 



(2)- 



99734647126127376 7 



459566406257( 2 ) - 249336617815318447' 
748009853445955327 



6565234375 7 ( 4 ) 



-Poisson 



-Poisson 
#2,8 



-Poisson 

yi,8 



-Poisson 

yo,8 



^Poisson 



2992039413783821287 



149601970689191064 7 - 229783203125 7(2) + 65652343750 7 ( 4 ) ' 

69&J 

5152958606448872767- 781262890625 7(2) + 1838265625007( 4 ) ' 

5984078827567642567 

274270831076537784 7 - 413609765625 7(2) + 91913281250' 
37400492672297766 
45956640625 ' 
274270831076537784 7 - 4136097656257(2) + 91913281250 7 ( 4 ) 
919132812507 ' 



while in all other cases where j < k, we have 

Appendix E: Relation to the blocking 
(disconnected) and connected correlation 
functions of fluids in disordered porous media 

In the statistical mechanics of fluids in a disordered 
porous medium, and more generally in the theory of 
disordered systems, the existence of two different types 
of averages, a "quenched" average over disorder and an 
"annealed" average over the equilibrated fluid, requires 
the introduction of several distinct correlation functions. 
At the pair level, these correlation functions are known 
as "connected" and "blocking" 72 i 73 (or "disconnected"). 
For a fluid they are defined according to 

P W(\t - r'|) = (p(r)p(r>)) - <p(r)> {p(r>)) (El) 

p[ 2) (|r-r'|) = (p(r)}(p(r')> ( E2 ) 

where p(r) is the microscopic fluid density within the 
open volume left by the disordered environment and the 



brackets denote an equilibrium in the presence of the 
latter. 

In the context of partially pinned fluid configurations, 
the same definitions apply provided that one associates 
the average in the presence of disorder with the condi- 
tional average ( ) e , the microscopic fluid density with 
p(r) (and thus the mean fluid density with (1 —c)p), and 
the average over the disordered environment with the 
double average over the pinned set and the equilibrium 
reference configuration. For instance, the blocking pair 
density function then reads 

p[ 2) (|r-r'|) = (<#r)) B <#r')) B ) . (E3) 

It is now easy to prove by following the line of reasoning 
as for deriving Eq. (|19p (see also Ref. 1 701 ) that 

(<tfr)) B (p(r')> B ) = (W)WYh) , (E4) 

so that the overlap Q c (oo) given in Eq. (fTB")) can also be 
written as 

Qc(oo) = 

{l _\ ?N j d *j dv'Q{a-\v-v'\)pf\\v-v'\). 

(E5) 
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After introducing the blocking total correlation func- 
tion through p^\r) = (1 — c) 2 p 2 [hb(r) + 1], one finally 
obtains that 



; (oo) - Qo(oo) = p dr Q(a - r)h b (r) 



(E6) 



The quantity that is relevant for extracting a point-to- 
set correlation length associated with the slowdown of 
dynamics is therefore simply related to the integral over 
a sphere of radius a of the blocking total correlation func- 
tion. (Note that even for hard spheres the latter does not 
equal —1 inside the core.) One could then envisage to use 
the approximations from liquid-state theory formulated 
in the context of fluids in disordered porous media^r— to 
compute the blocking function, although one may wonder 
if simple approximations, e.g., HNC or Percus-Yevick 66 , 
would be able to capture the behavior of the overlap. 
In addition, it should be recalled that whereas the solid 
porous matrix that leads to the average over the dis- 
order does not change with the thermodynamic point, 
the disorder associated with a pinned sets of particles in 
reference equilibrium fluid configurations does vary with 
packing fraction (and temperature in general). 
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